Dynamic properties of the spin-1/2 XY chain with three-site 

interactions 



Taras Krokhmalskii 1 ' 2 , Oleg Derzhko 1,2 , Joachim Stolze 2 , and Taras Verkholyak 1 ' 2 
1 Institute for Condensed Matter Physics, 
National Academy of Sciences of Ukraine, 
1 Svientsitskii Street, L'viv-11, 79011, Ukraine 
2 Institut fur Physik, Technische Universitat Dortmund, ^221, Dortmund, Germany 

(Dated: December 25, 2008) 

Abstract 

We consider a spin-1/2 XY chain in a transverse (z) field with multi-site interactions. The 
additional terms introduced into the Hamiltonian involve products of spin components related to 
three adjacent sites. A Jordan- Wigner transformation leads to a simple bilinear Fermi form for 
the resulting Hamiltonian and hence the spin model admits a rigorous analysis. We point out 
the close relationships between several variants of the model which were discussed separately in 
previous studies. The ground-state phases (ferromagnet and two kinds of spin liquid) of the model 
are reflected in the dynamic structure factors of the spin chains, which are the main focus in this 
study. First we consider the zz dynamic structure factor reporting for this quantity a closed-form 
expression and analyzing the properties of the two-fermion (particle-hole) excitation continuum 
which governs the dynamics of transverse spin component fluctuations and of some other local 
operator fluctuations. Then we examine the xx dynamic structure factor which is governed by 
many-fermion excitations, reporting both analytical and numerical results. We discuss some easily 
recognized features of the dynamic structure factors which are signatures for the presence of the 
three-site interactions. 

PACS numbers: 75.10.Jm; 75.40. Gb 
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I. MULTI-SITE INTERACTIONS AND JORDAN- WIGNER FERMIONIZATION 



Spin- 1/2 XY chains provide an excellent ground for studying various properties of quan- 
tum many-particle systems since after performing the Jordan- Wigner transformation these 
spin models can be reduced to systems of noninteracting spinless fermions [1]. One inter- 
esting issue which has emerged recently in the theory of quantum spin systems is the study 
of effects of multi-site interspin interactions. Such interactions may arise, e.g., in an effec- 
tive spin model for the standard Hubbard model at half filling in higher orders (beyond 
t 2 /U) of the strong-coupling t/U expansion J2j. Another example is provided by quan- 
tum spin systems with energy currents 0, 0, [J. As early as 1971 [6] M. Suzuki proposed 
generalized one- dimensional XY models with multi-site interactions, allowing for rigorous 
analysis by the Jordan- Wigner fermionization approach. An exactly integrable spin- 1/2 
XXZ quantum spin chain with three-site interactions was suggested in Ref. Q (see also 
Ref. [s). An XY version of this model was considered independently in Ref. |^. An- 
other version of the spin- 1/2 XY chain with three-site interactions was suggested in Refs. 
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Later on spin-1/2 XY chains with three-site interactions were considered in a num- 
ber of pap ers concerning quantum phase transitions, transport properties and entanglement 
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20|. Recently, in Ref. [2l| the spin-1/2 X Y chain with alternat 



ing three-site interaction has been introduced, whereas in Ref. [22| dynamic characteristics 
of a few quantum spin chains with multi-site interactions have been discussed. However, an 
exhaustive study of the dynamic properties of spin-1/2 XY chains with multi-site interac- 



tions, similar to that for conventional XY chains 
With our paper we attempt to fill this gap. 
In what follows we consider the Hamiltonian 
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251 ]. has not been performed yet. 



H — (j n (s n s n+1 + s v n s y n+A ) + D n (s„s 



y 

n+l 



y a x 



) 



+K n (s n S n+1 S n+2 + S n S n+l S n+2) + { S n S n+l S n+2 S n S n+l S n+2) 



Here J n and D n are the isotropic XY (or XX) exchange interaction and the ^-component of 
the Dzyaloshinskii-Moriya interaction between the neighboring sites n and n+l, respectively. 



K„ anc 



in Ref. 



E n are two types of three-site exchange interactions introduced in Refs. 



9 , 



respectively (see also Ref. 
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and 



26l where the general Hamiltonian (11. ip was introduced 



show the close relations between the models of Ref. (i 
three-site interactions) and of Refs. 
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as well). Q n is the transverse (z) external magnetic field at the site n. The sum in Eq. 
(11. ip runs over all N lattice sites; boundary conditions (open or periodic) are not important 
for the quantities to be calculated in the thermodynamic limit N — > oo. In this study in 
most cases we restrict ourselves to homogeneous chains with site-independent values of the 
interspin interaction constants and field, i.e. J n = J etc. 

We start by discussing the symmetry properties of the Hamiltonian (11.11) in order to 

s of Ref. |^ (i.e. with the XZY - YZX type of 
i.e. with the XZX + Y ZY type of three-site 
interactions) which were not discussed before. Consider a local spin rotation around the z 
axis 

S n K = S n COS K + S V n Sin <f) n , 
S l -> S n = ~ S n Sin <Pn + S V n COS (fi n , 

s n —> s n = s n . (1-2) 

Under that transformation, the parameters for the interspin interactions in the Hamiltonian 
(II .ip are mapped as follows: 

Jn -> Jn = Jn COS ((fi n +l - n ) + D n Sm (0 n+1 - n ) , 

D n -> D n = - J n sin (4> n+1 - 4> n ) + D n cos (<p n+1 - <p n ) , 
K n -> K n = K n cos (4> n+2 - 4> n ) + E n sin (<j) n+2 - <p n ) , 
E n -> E n = -K n sin (0 n+2 - <f> n ) + E n cos (0 n+2 - <p n ) . (1.3) 

This shows clearly that the rotations (II. 2p may be employed to simplify the Hamiltonian 
(II. ip by eliminating some of the interactions. For example, in a homogeneous chain we 



may achieve D = by setting 0„+i — <p n = ip with tan<^ = D/J 27|. The remaining 
coupling constants then are J = sgn(J) y/ J 2 + D 2 , K = ((J 2 — D 2 ) K + 2JDE) / (J 2 + L> 2 ), 
E = (-2JDK + (J 2 - D 2 ) / (J 2 + D 2 ). 

More interestingly, using Eqs. (II. 2p . (11.31) we can eliminate from the Hamiltonian Eq. 
( 11 .ip either of the three-site interactions, K n or E n . Introducing 6 n = (f) n +2 — (fin we note 
that K n = if tan^ n = —K n /E n . These 9 n values can be used for calculating the J n and 
D n by using <fi 2m - <p2m-i = <fi>2 - <fii + YJj=i - Qij-i) and <fi 2m +i - <fi>2m = 02m-i ~ 
02m + 02m-i; the surviving three-site coupling is E n = sgn(E n ) 'E 2 + K 2 . If, on the other 
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hand, we put tan6 n = E n /K n we get E n = whereas K n = sgn(K n ) a/ K% + E 2 . For the 
uniform chain, on which we mainly focus in what follows, 9 n = 8, 2m — 02m- 1 — 02 ~ <fii, 
02m+i — 02m = # — 2 + 0i, and assuming 2 — 0i = 9/2 we get J = J cos(#/2) + D sin(0/2), 
D = -Jsin(0/2) + Dcos(6/2) and either K = 0, E = sgn(E)VK 2 + E 2 if tan# = —K/E 



ot E = 0, K = sgn(K)VK 2 + E 2 if tan# = E/K. 

To summarize this part, we have shown that while studying effects of three-site interac- 
tions on the basis of the model (11.11) it would be sufficient to consider the model ( 11. ip with 
either nonzero parameters J n , D n , K n or J n , D n , E n since the most general case when all 
four types of interactions have nonzero values can be reduced either to the former case or 
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can be 



to the latter case. One consequence of this is that the model considered in Ref. 
reduced to the model considered in Ref. Namely, starting from the model with J ^ 0, 
K ^ 0, D = E = and choosing n + 2 — <Pn = vr/2, n +i — 4> n = vr/4, n = nir/A we arrive 
at the model with J = J/y/2, D = -J/V2, E = -K, K = 0. Vice versa, starting from the 
model with J^O, D ^ 0, E ^ 0, K = and performing the same transformation we arrive 
at the model with J = (J + D) /y/2, D = {—J + D) /y/2, K = E, E = 0. In our study of 
the dynamic properties of quantum spin chains with three-site interactions we focus on the 
case J ^ 0, K ^ 0, D = E = leaving a detailed study of other cases for the future. 

The peculiar nature of the three-site interactions studied here becomes clear after per- 
forming the Jordan- Wigner transformation 

s + =s x + is^ = P ir t s - = s x _ ; s « = p , r 

C n Pn— l^n j C n P n —\S n , 
m m 

p m = n (i - acjc,-) = n (-24) • a-*) 

Plugging (11.41) into (ll.ip we find 



tt _ ( J n + iDn t 1 Jn-lDn J 

n — 2^ I — 2 — 2" — 



. C n C n+2 . C n+2 C n 



" n »( c ic»-^)). (1-5) 



i.e. the Hamiltonian of the spin model is a simple bilinear form in terms of spinless fermions. 
For the uniform case it is convenient to employ periodic boundary conditions in the spin 
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Hamiltonian (11.11) . That leads either to periodic or antiperiodic boundary conditions in the 
fermion Hamiltonian (ll.5p . depending on whether the number of spinless fermions is even 
or odd. Either the nearest-neighbor hopping integrals or the next-nearest-neighbor hopping 
integrals can be made real (or purely imaginary) by applying a gauge transformation, which 
is the analog of Eqs. (jl.2p . (11.31) in the spinless- fermion picture. 

In the uniform case we can diagonalize (11.51) by performing the Fourier transformation 



7= exp c «' ° n = ~7T? exp (- [Kn ) c ^ I 1 - 6 ) 



with k = (2n/N)m [k = (2n/N)(m + 1/2)] in the subspaces with odd [even] numbers of 
spinless fermions and m = —N/2, —N/2 + 1, . . . , N/2 — 1 (if N is even) or m = —(N — 
l)/2, —(N - l)/2 + 1, . . . , (N - l)/2 (if N is odd) arriving at 



K E 

A K = J cos k + D sin k — — cos(2k) — — sin(2/t) + Q. (1-7) 

From Eq. (11.71) we immediately conclude that the external magnetic field Q plays the role 
of a chemical potential for spinless fermions. More interestingly, the two terms proportional 
to cos(2/t) and sin(2/t) in the elementary excitation spectrum A K (1 1.7ft which arise from the 
three-site spin interactions, may modify A K drastically, leading to new ground-state phases. 
We consider the case J ^ 0, K ^ 0, D = E = 0; for this case A_ K = A re . As long as 
| if | < | J |/2 the spinless fermion system may possess only two Fermi points, whereas for 
sufficiently strong three-site interaction, \K\ > \ J\/2, the spinless fermion system may also 
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is shown 



possess four Fermi points. The ground-state phase diagram obtained in Ref. 
in Fig. [U Two different spin liquid phases reflect the importance of the existence of two 



versus four Fermi points, as discussed already earlier [TjJ. As a result, in addition to the 
conventional quantum phase transition between the spin liquid I phase and the ferromagnetic 
phase, the spin model may exhibit also quantum phase transitions 1) between the spin 
liquid I phase and the spin liquid II phase and 2) between the spin liquid II phase and the 
ferromagnetic phase, as well as a point where three phases meet (see Fig. [p. The spin liquid 
I phase and the spin liquid II phase are characterized by a change in the pow er-law decay 
and oscillating factor of spin-spin correlations (see Eqs. (20) - (23) in Ref. Ill 



To describe a zero-field quantum phase transition from the spin liquid I to the spin liquid 



II phase at \K\ = 2\J\, I. Titvinidze and G. I. Japaridze 



111 ] introduced the order parameter 



-2.5 
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FIG. 1: The ground-state phase diag 



D = E = discussed earlier in Ref . 



ram of the homogeneous model (jl.ip with J = ±1, K ^ 0, 
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J. The region — l+i^/2 < < l+if/2 corresponds to the spin 
liquid I phase (two Fermi points) [dark-gray], the regions K < —1/2, l+K/2 < < — K/2— 1/(4K) 
and -RT > 1/2, —K/2 — 1/(4K) < < — 1 + K/2 correspond to the spin liquid II phase (four Fermi 
points) [light-gray], the remaining regions correspond to the ferromagnetic phase [light]. The black 
dots correspond to the sets of parameters which we most often use below to discuss dynamic 
quantities. 



7] constructed from the average length £ of the ferromagnetic string in the ground state 
(see Eqs. (31) - (33) in Ref. Illl ). On the other hand, P. Lou et al. in their study 
of the (homogeneous) model (11.11) with J^O, E ^ 0, D = K = introduced the scalar 
chirality parameter O = — (1/iV) J2 n ( s n-i s n s n+i ~ s n-i s n s n+i) an d calculated this quantity 
in the ground state at Q = (see Eqs. (10), (11) in Ref. Il2l ). A nonzero value of O 
signalizes the appearance of a new (chiral) spin liquid phase which emerges when three-site 
interactions exceed a critical value. Eliminating XYZ — YZX terms from the Hamiltonian 
by the unitary transformation discussed above one arrives at the spin model with J = J / a/2, 
D = — J/y/2, K = E, E = (i.e. the model similar to the one considered in Ref. Illl ) and 
0^0 = —(1/N) J2 n ( s n-i s n s n+i + s n-i s n 3 n+^L The latter quantity O apparently has no 
relation to the order parameter rj used in Ref. Ill 



It is also worth mentioning here that the spin model (11 .ip can be written as a one- 
dimensional model of hard-core bosons after introducing the on-site creation and annihilation 



operators s+ = + is^ and s~ 



a® i qV 



1 S y S 
la ra> °n 



S n S n 



1/2. The hard-core boson model 



is obtained by taking the U — > oo limit of the boson Hubbard model. With this mapping 
in the case of a conventional transverse XX chain the ferromagnetic phase corresponds to 
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the Mott insulator with (1/N) J2 n ( s n s n) = or (1/iV) J2n( s n s n) = 1 (^ ne thermodynamic 
average taken at zero temperature, T = 0), whereas the spin liquid (spin liquid I) phase is 
the superfluid with < (1/N) J2 n ( s n s n) < 1- As a function of the field / chemical potential 
Q the model displays two superfluid to Mott insulator transitions at Q = ±| J\ |28|]. After 
switching on the three-site interaction K ^ 0, the picture remains qualitatively the same as 
long as \K/J\ < 1/2. If \K/J\ exceeds 1/2 we face an additional transition which manifests 
itself as an extra cusp in the dependence of (1/N) J2 n ( s n s n) on ^- Thus, following the 
ground-state average boson number per site (or m z = (1/N) J2 n ( s n) = En(4 s ») — 

1/2) as a function of Q one may reproduce the various phases and the phase transitions 
between them shown in Fig. [TJ An alternative way to follow the changes in the ground-state 
dependence of m z on Q is to examine the ground-state susceptibility Xzz = dm z /dQ as a 
function of Q. We notice that the ground-state dependence of — Xzz on Q is the same as that 
of p(E = 0) on tt, where p(E) = (1/N) J2 K S(E-A K ), A K = J cos k-(K/2) cos(2/t)+fi is the 
one-particle density of states. As a result, the ground-state dependence Xzz v s exhibits a 
square-root van Hove singularity along the lines separating different phases in Fig. [TJ (The 
only exception are the two points K = ±| J\/2, Q = =f|J| + K/2 at which A K cx k 4 and 
therefore Xzz displays a van Hove singularity with the exponent 3/4.) The divergence of the 
uniform static zz susceptibility implies a "ferromagnetic" character of the associated phase 
transitions. 

To summarize, there is no doubt that while "the two Fermi points spinless fermions" 
transform into "the four Fermi points spinless fermions" some noticeable changes in the 
properties of the spin model should take place, however, a transparent quantity associated 
with this modification of the Fermi surface topology which may play the role of the order 
parameter is still lacking. 

Finally, it is worth noting the studies on the one-dimensional Hubbard model with next- 
nearest-neighbor hopping because the noninteracting limit of that model resembles Eq. (\1.7\i 
(see, e.g., Ref. [29| and references therein). In contrast to those studies, we calculate two- 
and many-particle correlation functions (although for a system of noninteracting spinless 
fermions (11.71) ) which are related to two-spin correlation functions for an interacting quantum 
spin system. 

In our study of dynamic properties of the spin model we focus on the dynamic structure 
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factor 

N ,00 

S AB (K,cu) = J2 e xv(-™ 1 ) / dtexp (kot)((A n (t) -(A)) (B n+l (0) -(B))), (1.8) 

where A n , B n are some local operators attached to the site n (like s", a = 
x,y,z or = s x n s^ +1 + s y n s y n+1 etc.), A n (t) = exp(iift)A n exp(-iHt), ((...)) = 

Tr(exp(— [3H){. . -))/Tr exp(— /3H), (A) = (l/N)^2 n (A n ). Knowing the dynamic structure 
factors we can find the corresponding dynamic susceptibilities according to well known re- 
lations (see, e.g., Ref. |30T ) . 

In the next section (Sec. [TT]) we report a closed-form expression for the zz dynamic 
structure factor S zz (k,u>) (i.e. A n = B n = s z n in (11.81) ) and for some similar dynamic 
structure factors Si(k,u) (see Eqs. (12. 3p . (12. 4p below) all of which are governed by two- 
fermion (particle-hole) excitations. Then, in Sec. IHH we discuss general properties of the 
two-fermion excitation continuum focusing on spectral boundaries, soft modes, singularities 
etc. We also contrast generic and specific features of various two-fermion dynamic quantities. 
In Sec. IIVI we examine many-fermion dynamic quantities focusing, in particular, on the xx 
dynamic structure factor S xx (k,u) (i.e. A n = B n = in (11.81) ). We report exact analytical 
results 1) in the high-temperature regime (3 = (T — > 00) and 2) in the strong-field regime 
in the ground state as well as precise numerical results for arbitrary temperatures. Finally, 
in Sec. |V] we summarize our findings. Some selected results of the present study were 
announced in a conference paper [31 ]. 



II. zz DYNAMIC STRUCTURE FACTOR AND SOME OTHER TWO-FERMION 
DYNAMIC STRUCTURE FACTORS 



We start with the calculation of the transverse dynamic structure factor S zz (k,u) that 
corresponds to A n = B n = in Eq. (11.81) . Since according to Eq. (11.41) sj — c^c n — 1/2 
the calculation of the transverse dynamic structure factor is very simple 23j, |32|. Using Eqs. 
( II. 6p . (I1.7P and the Wick-Bloch-de Dominicis theorem we end up with the result 

/TT 
dK!Ti Kl (1 - n Kl+K ) 8(u> + A K1 - A Kl+K ) 
TT 

n Kl (1 — n Kl+K ) 



E 



—f _2_ (A - A 



(2.1) 



Here n K = 1/ (1 + exp(/3A K )) is the Fermi function, and {k*} are the solutions of the equation 



it) + A K * - A K * +K = 0. 



(2.2) 



There are more local spin operators which in fermionic representation are given by a 
product of two Fermi operators, 

^(l) — o x o x _l ~y q v — 1 fj r , _ r J \ 



x Q y 

n°n+l 



S n S n+l 2 (^ C n C n+l + c n c n 
M) — „x z x , y z V _ _}_ ( \ _ t 



+1 ' 



+(2) _ x z y _ y z x __Af r t r i r 
b n ~ a n a n+l^n+2 °n n+l°n+2 ~ ^ \^n u n+2 i u n<*- n + 2 



,t 



(2.3) 



etc. The correlation functions (d { n ] (t)d%(0)) , (d { n ] {t)d%(0)) , (^(^^(O)), 

(2) (2) 

(tn {t)tn+i (0)) an d the corresponding structure factors Sj(k,u)), Sd(k,u)), a>), Se(k,u)) 

therefore can be again easily calculated with the result 



d.KiBj(K,i, Ki + k)C(k,i, Ki + k)S (uj — E(K\, Ki + k)) 

Bj(ki, Ki + k)C(ki, Ki + k) 



8k 



-E(kx, ki + k) 



2 Ki + K 2 



Bj(ki, k 2 ) = COS 

r> / \ . 2 K l + K 2 

B d (ki,k 2 ) = sin — , 

B k (ki, K 2 ) = | cos 2 («i + K 2 ) , 

5 b (aci, « 2 ) = i sin 2 (ki + k 2 ) , 
C(«i, « 2 ) = n Kl (1 - n K2 ) , 
E(ki, k 2 ) = -A Kl + A K2 . 



(2.4) 



We note that S zz (k,uj) (12. ip is also given by Eq. (12.41) with B zz (ki, k 2 ) = 1. 

In Figs. [2J [31 IH we show the gray-scale plots for the different two-fermion dynamic 
structure factors (12.11) . (12.4p for several representative sets of the Hamiltonian parameters, 
at zero and infinite temperatures. Comparing the closed-form expressions (12.11) . ( 12. 4ft and 
the gray-scale plots in Figs. [21 [31 0] we conclude that 1) the generic properties of all these 
two-fermion dynamic quantities are controlled by the 5-function containing E(ki, k 2 ) in Eqs. 



FIG. 2: S zz (k,u) for the model {H) with J = 1, D = E = 0, K = 0.5, n = (panel a), K = 2, 
Q = (panel b), # = 2.5, O = (panel c), K = 0.5, O = 1.25 (panel d), K = 2, O = -1.125 
(panel e), if = 0.5, f2 = —0.75 (panel f). Panels a, b, c refer to the ground state (T = 0), panels 
d, e, f refer to the low temperature [3 = 10. 

(12. ip . (I2.4p in the high-temperature limit T — ► oo (/3 — > 0), whereas in the low-temperature 
limit T — > (/3 — > oo) the factor containing Fermi functions, C(k\, k 2 ), becomes important 
in addition and 2) the specific properties of the various two-fermion dynamic quantities 
arise only owing to different functions Bj(ki, k 2 ). In the next section we further examine 
the two-fermion dynamic structure factors, revealing their similarities and contrasting their 
differences. 



III. TWO-FERMION EXCITATION CONTINUUM: GENERIC VERSUS SPE- 
CIFIC PROPERTIES 

Let us discuss the properties of the two-fermion excitation continuum which is probed by 
a number of dynamic quantities like the transverse dynamic structure factor S zz (k,u), the 
dimer dynamic structure factor Sj(k,uj) etc. For the model under consideration here, (II. ip 
with coupling constants D = E = 0, the elementary excitation spectrum (11.71) differs from 
that of the standard homogeneous XX chain by the cos(2fi;)-term. That term has important 



23 



consequences which we explore by generalizing the work of J. H. Taylor and G. Muller 
on conventional XY chains. 

We start with the high-temperature limit T — >• oo. A two-fermion dynamic quantity (I2.4p 
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FIG. 3: Two-fermion dynamic structure factors Sj(k,uj) (panel a), Sd(k,uj) (panel b), Sk(k,w) 
(panel c), S e {k, u) (panel d) for the model dTTTJ with J = 1, D = E = 0, K = 2.5, O = at T = 0. 




FIG. 4: Two-fermion dynamic structure factors S zz (k,uj) (panel a), Sj(k,uj) (panel b), Sk(k,uj) 
(panel c) for the model (jl.ip with J = 1, D = E = 0, K = 2.5 at T ^ oo. Note that for infinite 
temperature the two-fermion dynamic structure factors are field- independent. 
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may have a nonzero value at a point (k, uj) in the wavevector - frequency plane (we assume 

UJ > 0, — 7T < K < 71") if 

cj = £■(«!, k 2 ) = -A Kl + A K2 , k = -Hi + K 2 (mod(27r)), 

K 

A K = Jcosk cos(2k) + Vt, (3.1) 

2 

where — tt < K\ < tt. We rewrite the function E(k\, k\ + k) in the form 

E(k 1 ,k 1 + k) = 2 sin ^ sin ^ + k^ ( _ ^ + 2if cos^cos ^ + k^ (3.2) 

and solve the equation 8E{ki,ki + k)/8ki = with respect to K\, or more precisely with 
respect to x = cos (k/2 + Ki), to find 



x± = ttftzztt; ± +1 (3-3) 



8iTcosf y \8Kcos^J 2 

For < 1/2 there is only one pair of Ki which solves (13.31) and fulfills the condition 

|x| < 1, thus yielding a stationary point of the function E(k±, Ki + k). That pair is given by 
k± = arccosx - — k/2, k^ = — arccosx - — k/2 for JK > 0, and by kf = arccosx + — k/2, 
k~{ = — arccosa;" 1 " — k/2 for JK < 0. In the opposite case \K/J\ > 1/2 there are two 
such pairs, k\ = arccosx - — k/2, k± = — arccosx" — k/2 and kf = arccosx + — k/2, 
k± = — arccosx + — k/2. As a result, the upper boundary oj u (k) of the two-fermion excitation 
continuum is given by 

= \E(ki,ki +k)\, JK>0 (3.4) 

or 

Q+(k) = \E(k+,kf + k)\, JK<0. (3.5) 

Note that for \K/ J\ > 1/2 the first derivative of E(ki, ki + k) with respect to K\ is zero also 
along Q + (k) < Q-(k) = uj u (k) for JK > or along f2_(/t) < Q + (k) = uj u (k) for JK < 0. 

Due to the presence of a two-particle density of states the two-fermion dynamic quantities 
(12. ip and (12 .4p may exhibit a van Hove singularity along the lines 

W a («)={fl_(«),fi + («)}, (3.6) 

or more precisely along Q^(k) [Q + (k)} if JK > [JK < 0] for \K/J\ < 1/2 and along both 
lines f2_(/t) and Q + (k) for \K/J\ > 1/2. Thus, a sufficiently strong three-site interaction 
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FIG. 5: Different kinds of potential van Hove singularities of the two-fermion dynamic quantities. 
Frequency profiles of S zz (k, uj) (i.e. Si(k,uj) with Bi(ki,K2) = 1) at T — ► oo, J = 1. Panel a: 
K = 7r/2, X = 0.25 (dotted line), K = 0.5 (dashed line), K = 2.5 (solid line). Panels b and c: 
K = 2.5, k = 2.70 (dotted lines), n = k° 2.73887681 (dashed lines), k = 2.78 (solid lines). The 
initial slope of the dashed line in panel c corresponds to the dependence oc u~ 2 / 3 . 

K increases the number of van Hove singularities. This is nicely seen in Fig. H] (and 
also in Fig. [5]) where the two-fermion dynamic structure factors for J = 1, K = 2.5 at 
T — > oo are plotted. [Recall that B zz (k±, K2) = 1 and therefore Fig. with gray-scale 
plot for S zz (k,u) most transparently demonstrates a new (low-frequency) line of van Hove 
singularities emerging for \K/J\ > 1/2.] 

Interestingly, in addition to the conventional inverse square-root van Hove singular- 
ity a singularity with exponent —2/3 may occur as \K/J\ > 1/2. In fact, plotting 
d?E(Ki, K\ + /c)/c?«i| R =ft vs k, we note that 3 2 E(ki, k,\ + K)/dKf\ Ki=ki = for k = k° 
which satisfies x + = 1, JK > (or x~ = —1, JK < 0). A similar analysis of 
d 3 E(ni, Ki + K)/dn\ \ k vs k shows that d 3 E(ni, Ki + k)/9k 3 | Ki= - i 7^ for k = k°. More- 
over, we find that E(k\, K\ + k) = for these values of k and K\. This immediately implies 
that Si(K ,e) oc e~ 2//3 , e — > +0. In Fig. Owe demonstrate potential van Hove singulari- 
ties of the two-fermion dynamic structure factors. In particular, we illustrate the van Hove 
singularity with exponent —2/3 for J = 1, K = 2.5. Solving equation x + = 1 (13. 3p with 
respect to k we find /t w 2.73887681. The frequency profiles around this value of wavevector 
clearly show two types of van Hove singularity, i.e. with exponent —1/2 (in most cases when 
uj — > uj s (k) — e; e.g., the dotted and solid lines in Figs. [5b and[5b) and with exponent —2/3 
(only when n = and uj —>■ e; the dashed lines in Figs. [5b and[5t). 

We consider now the case of zero temperature T = when the Fermi functions entering 
the function C(«i,«2) in Eq. (12.41) become extremely important. They imply that in the 



13 



ground state we have to require in addition to Eq. (13. ip A Kl < and A K2 > 0. We first 
consider the case \Kf J\ < 1/2. Plotting the dependence of E(k\, ki + k)C(ki, ki + n) on Ki 
we find that the two characteristic curves, A Kl+K with K\ satisfying A K1 = (or E(k\, k% + k) 
with k>i satisfying A Kl = 0) and — A Kl with K\ satisfying A K1+K = (or E(ki,Ki + k) with 
K\ satisfying A Kl+K = 0), play a special role. Solving the equation A^ = for k, or more 
precisely for y = cos k, we find 

y ± = & ± S+FI- < 3 - 7 > 

Taking into account that \y\ < 1, we see that for \K/J\ < 1/2, Eq. (13. 7p may yield two k 
values, k = arccos y~ , k = — arccos y~ (for JK > 0), or k = arccos y + , k = — arccos y + (for 
JK < 0). For \Kj J\ > 1/2 there may be two such pairs of k, k = arccos y~, k = — arccos 
and k = arccos y + , k = — arccos y + ; we will discuss that case later. 
For \K/J\ < 1/2 we consider two characteristic lines 

U)_{k) = |A al . ccosy - +K |, UJ_(k) = |A_ arccos y-+ K I) JK > (3.8) 

or 

= l^-arccosy+4-Kl) UJ +{ K ') = I ^- arccos y++n, \ j JK < 0. (3-9) 

The smaller one of the two values u 1 _(k), i = — , + for JK > [uj 1 + (k), i = — , + 
for JK < 0] gives the lower boundary of the ground-state two-fermion excitation con- 
tinuum ui(k), whereas the other (larger) one gives either the upper boundary of the 
ground-state two-fermion excitation continuum u u (k) or the middle boundary of the 
ground-state two-fermion excitation continuum u m (K). The former case occurs if k\ 
yielding dE{K%, k% + k) / dn\\ - = (see above) belongs to the region of k± where 
E(ki,Ki + k)C , (k 1 ,k 1 + k) = (as, e.g., is seen in Fig. [6b for small \k\). In the lat- 
ter case, when k x yielding 8E(ki,ki + K)/dK\\ ~ = does belong to the region of Ki 
where E(k\,k\ + «)C(«i,«i + k) ^ 0, the larger value of u % _{k), JK > (13.81) [u)+(k), 
JK < (13. 9p ] gives the middle boundary of the ground-state two-fermion excitation con- 
tinuum uj m {K) whereas the upper boundary uj u (k) is given by Eq. (13.41) [Eq. (13.51) ]. For the 
frequencies to between uji(k) and u m (K) [uj m (K) and uj u (k)] equation (I2.2p has one solution 
[two solutions] k\ . Thus, the ground-state Si(k,u>) changes by a factor 2 at the middle 
boundary u) m (K). 
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The soft modes kq can be determined from the equation u^(kq) = 0, JK > [k>+(«o) — 0, 
JK < 0]. Therefore, if Jif > 0, 

k = {0, ±2 arccos y~ } , y~ > 0, 

Ko = {0,=f2 arccos y~ ±2vr} , y~ < (3.10) 

or, if JK < 0, 

kq = {0, ±2 arccos y + } , y + > 0, 

Ko = {0,^2 arccos y + ±27r}, y+ < 0. (3.11) 

We now turn to the case \K/J\ > 1/2. As already mentioned, equation (13.71) may 
yield two pairs of k, k — arccos?/ - , k = — arccos y~ and k = arccos y + , k = — arccos y + , 
and therefore all four characteristic lines in the k-u plane, w_(k), ujZ(k), u+(k,), uj+(k) 
given by Eqs. (13.81) . (I3.9P come into play simultaneously. Thus, as \K\ exceeds \ J\/2 an 
"extra" ground-state two-fermion excitation continuum emerges. Its lower boundary and 
upper /middle boundary are given by formulas (13.81) . (13.91) [in the case when Eqs. (13.81) . (13. 9p 
give the middle boundary, the upper boundary is given by one of the formulas in Eqs. (13.41) . 
(13. 5ft ]. The number of soft modes increases but cannot exceed 9. 

In the vicinity of a soft mode kq the lower boundary of the two-fermion continuum in most 
cases displays a "V" shape, i.e. it is proportional to \k — Kq\. However, it is worth noting that 
a parabolic shape oc (k — k>o) 2 is also possible for suitable parameter combinations. To see 
this we recall that the lower boundary is basically determined by the equation = with 
Afe = J cosk — {K/2) cos(2/c) + Q. For the conventional XX chain (K = 0) in a transverse 
field Q we see that at the critical field values Q = ±|J| the dispersion A& oc (k — ko) 2 , 
with a corresponding parabolic shape of the lower continuum boundary near a soft mode at 



kq. [This can be seen, e.g., in Fig. llf of Ref. 24|.] At the critical field the nature of the 
ground state of the XX chain changes from partially filled to completely filled or completely 
empty in terms of Jordan- Wigner fermions. For nonzero K the cos(2/c)-term may generate 
additional maxima or minima in Afc. That implies the emergence of additional critical 
field values, at which k regions near the additional minima or maxima of A& open or close 
for occupation by Jordan- Wigner fermions. These critical values correspond to the lines 
separating different ground-state phases in Fig. [TJ Along these lines we expect parabolic 
behavior of the lower two-fermion continuum boundary (see, e.g., Figs. [2b, [2H, [2b, [2f). 
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FIG. 6: Characteristic lines of the two-fermion excitation continuum. We assume J = 1 and 
K = 0.5, = 0, = -0.5, = 1 (panels a, b, c), K = 2, O = 0, O = -0.5, = 1 (panels d, 
e, f), IT = 2.5, O = 0, = —0.5, = 1 (panels g, h, i). x and + symbols correspond to 0_(/c) 
and + (k), respectively. Solid, dashed, dash-dotted and dotted lines correspond to u~z(k), loZ(k), 
u}\{k) and ujZ(k), respectively. (The panels a, d, g can be compared to the panels a, b, c in Fig. 

El) 



Summarizing this part, we report in Fig. [B]all the characteristic lines of the two-fermion 
excitation continuum discussed above. In these plots symbols correspond to fi_(«;) (x sym- 
bols) and Q + (k) (+ symbols), whereas lines correspond to oj1(k) (solid), uZ(k) (dashed), 
lo\{k) (dash-dotted), uj+(k) (dotted). For \K/J\ < 1/2 only three characteristic lines are 
relevant, but if \K/J\ > 1/2 all six lines are relevant. These lines are important not only for 
understanding the distribution of the two-fermion dynamic structure factors Si(k,u) ( 12. II) . 
( 12.41 over the k - u plane (see gray-scale plots in Figs. El El SJ) but also for many-fermion 
dynamic structure factors (like S xx (k,u>)) at low temperatures as will be discussed in the 
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FIG. 7: Illustrating the role of I?-functions. S zz (tt/2,lo) (panel a), Sj(ir/2,ui) (panel b), Sk(^/2,u>) 
(panel c) at T — > oo. J = 1, K = (dotted lines), iT = 0.5 (dashed lines), K = 1 (solid lines). 

next section (see gray-scale plots in Fig. HD1) . 

Finally, we discuss some specific properties of the two-fermion dynamic structure factors 
( 12. 4p controlled by different 5-functions. Comparing different panels in Figs. [3] and H] we 
observe a number of small but definite differences for the detailed distributions of Si(k,u) 
over the k — u plane. 

To be specific, we may focus on the dynamic dimer structure factor Sj(k,u). It is 
known that Sj(k, uj) does not diverge along the upper boundary owing to Bj(k x , k 2 ) for the 
conventional XX chain, i.e. when K = (see, e.g., Ref. |25|and references therein). This can 
be also seen in Fig. [7b where the dotted line corresponds to K = 0. This changes, however, 
if K ^ 0: the dynamic dimer structure factor exhibits a van Hove singularity along the 
upper boundary (the dashed and solid lines in Fig. [7_b) indicating the presence of nonzero 
three-site interactions. 

Next we may consider Sk(k,u). As can be seen in Figs. Hb and [7b, the van Hove 
singularity at the upper boundary becomes less distinctive as K increases. To explain 
this, we introduce the notation x = cos(k/2 + ki) and rewrite 8E(ki, k\ + K)/dKx ( 13. 2ft as 
2sin(/t/2)(2ircos(/t/2)(2x 2 - 1) - Jx), whereas B k {k x ,k x + k) (J22D as (l/4)(2x 2 - l) 2 . In 
the vicinity of the upper boundary the denominator in Eq. (12.41) for Sk{k, tends to zero, 
however, in the limit K — > oo the numerator in Eq. (12.41) for Sk{k, uj) becomes proportional 
to the denominator squared which makes the fraction equal to zero. Thus, for any finite 
large K the van Hove singularity at the upper boundary does exist (although with increasing 
K it is harder to find it numerically) and it disappears only in the limit K — ► oo. 
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IV. MANY-FERMION DYNAMIC QUANTITIES 



In this section we discuss many-fermion dynamic quantities fixing for concreteness our 
attention to the xx dynamic structure factor S xx (k,u). First we report two analytical re- 
sults referring to the high-temperature limit and to the zero-temperature strong-field regime, 
respectively, and then we turn to high precision numerical data for arbitrary values of tem- 
perature and the Hamiltonian parameters. 

We first consider the xx two-spin time-dependent correlation function (sj(t)sj +n ) at T — > 
oo. Since the Zeeman term commutes with the Hamiltonian of the considered model (11. ip 
and in the high-temperature limit exp(—j3H) — > 1 and consequently the averages of spin 
operators are zero we can easily extract the dependence on the transverse field fi 

(^m +n ) = cos (fit) (s*(t)s* +n )\ n=0 . (4.1) 

Thus we can proceed assuming fi = 0. Next we substitute into the spin correlation 
function on the r.h.s. in Eq. (14.11) a short-time expansion sj(t) = + i [H, sj] t — 



(1/2) [H, [H, s;-]] t 2 + . . . [33] and after simple but tedious calculations find 

4 * L-n = Ko 1 - + 3 o 1 )t 2 + ...), (4.2) 



J j+n/\n=o i \ 8 32 

where for generality we have considered a model with position-dependent couplings. Eq. 
(14. 2 p is consistent with the Gaussian decay 

Mm s i + n) = Snfl cos (fit) exp f- + A (4.3) 

for the model with position-independent couplings. Using MAPLE codes we checked that 
the terms up to t 4 in the short-time expansion for (sj(t)sj +n ) indeed agree with Eq. (14. 3p . 

Alternative (although not independent) arguments supporting Eq. (14. 3 p follow Refs. 
34], |35|. We examine the continued-fraction coefficients A& of the relaxation function 
c xx (z) = 4 f °°dtexp(-zt)(sj(t)sj) = l/(z+Ai/(z + A 2 /(z + . . .))) at T -> oo. The sequence 
of the continued-fraction coefficients reflects the time dependence of the associated au- 
tocorrelation function, and, in particular, when A& = kA then 4(sJ(t)s|) = exp(— At 2 /2). 
The sequence A& can be determined by the methods elaborated in Ref. [35!. We designed a 
MAPLE program which in a reasonable amount of time calculated A& for k = 1, 2, 3, 4 and 
confirmed the Gaussian decay (14.31) . 
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FIG. 8: Time dependence of the xx autocorrelation function in the high-temperature limit {(5 = 
0.0001); J = 1. Panel a: (s*(t)s*) at (1 = 0, K = 0, 0.5, 1, 1.5, 2, 2.5, 3 (from top to bottom). 
Panel b: {sj(t)Sj) at K = 1, O = 0, 0.5, 1, 1.5, 2 (from top to bottom). Symbols correspond to 
Eq. (|4.3|) ; lines correspond to numerical data obtained for N = 400, j = 41. 

Finally, from our calculations we also find a more general result for the homogeneous 
model (II. ip which is given by Eq. (0~2D after the substitution J 2 -> J 2 + D 2 , K 2 -> if 2 + £ 2 . 

To summarize, in Fig. [8] we compare analytical predictions according to Eq. (14.31) (sym- 
bols) with numerical calculations (lines) (see below) and observe an excellent agreement 
between both sets of data. Eq. (14.31) provides an extension of the well-known result for 
the conventional transverse XX chain [36] for the kind of three-site interactions considered 
here. The presented arguments in favor of the Gaussian decay (14.31) may be put even on a 
more rigorous foundation using the approach elaborated in Ref. |37. 



Next we turn to the zero-temperature strong-field regime. More precisely, we consider the 
ferromagnetic phase (light) in Fig. [TJ In the ferromagnetic phase the ground state of the spin 
model is completely polarized, i.e. |GS) = Yl n I i)n (|GS) = n n I T)n) f° r positive (negative) 
Q, that permits easily to take into account the Jordan- Wigner sign factors entering the 
formula for (s°!j(t)s* +n ) [38j. Assuming, for example, |GS) = Yl n I l)n (light region in the 
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upper half-plane in Fig. [T]) we immediately get 

K 

n^oo i r dreexp ( i ( Kn _ A(t t)) (4.4) 



with A K = Jcosk — (K/2) cos (2k) + > 0. For |GS) = n n I t)n (light region in the lower 
half-plane in Fig. [1]) we have 



4 = 4E ex P H «* =F tt) n - A„f)) 



AT-+oo 1 

^ 2^ 



AT 

r 

d^exp (— i ((« =f 7r) n - A K t)) (4.5) 



with A K = J cos k — (K/2) cos (2k) + fi < 0. Eqs. (|4~4|) . fT4~5|) contain the result for the 

= exp(-i|Q|0(-i)" UJt), where J„(z) 



conventional transverse XX chain [381], 4( s? 



is the Bessel function of the first kind 



39(. Some further properties of (sj(t)sj +n ) are 



collected in the Appendix. We only notice here that in the regime considered the xx time- 
dependent correlation function oscillates, with the envelope decaying proportional to t -1 / 2 
as t -> oo for \K/J\ ^ 1/2 or proportional to r 1/4 for |AT/J| = 1/2 (see Eqs. flA~6l) . flA~7l) . 
( TA~8|) and Fig. EJ). 

Eqs. (14 .3p and (14.41) . (14.51) immediately yield the xx dynamic structure factor (11.81) . In 
the high-temperature limit we have 



2 

i.e. the xx dynamic structure factor displays a ^-independent Gaussian ridge centered at 
frequency with the width controlled by interspin interactions. In the zero-temperature 
strong-field regime we have 

7T 

S xx (k, u) = -8 (u - lu («)) , 

+ Jcosk - yCOs(2«), Q > 0, 

|Q| + Jcosk + -jCos(2k), Q < 0, 

i.e. the xx dynamic structure factor displays a 5-peak along the line u*(k) (14.71) in the k 
- u) plane. Interestingly, if K ^ the symmetry of Eq. (14.71) with respect to the change 
Q — > — Q is broken, in agreement with the ground-state phase diagram shown in Fig. [TJ 
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FIG. 9: Time dependence of (sj(t)sj +n ) (panel a: n = 0; panel b: n = 1) in the low-temperature 
limit (/5 = 100); J = 1, K = 0.25, = 10. Lines correspond to numerical data for N = 400, 
j = 101. Symbols correspond to the long-time asymptotics given by Eq. (|A.6|) . (We notice that 
the asymptotic becomes accurate already for short times.) Solid (dashed) lines or filled (empty) 
symbols refer to real (imaginary) part of (sj(i)sj +n ). 

We turn next to the case of arbitrary values of temperature and the Hamiltonian param- 
eters. In this case we calculate the xx dynamic structure factor numerically. The numerical 
approach for calculating dynamic quantities was explained in detail earlier j^dl jjoT ] . To cal- 
culate (s|(t)s| +n ) we express the spin operators s x entering this quantity in terms of the 
Fermi operators c K , cj. according to (11.41) . (jl.6p obtaining as a result an average of a product 
of large number of Fermi operators attached not only to the sites j and j + n but to two 
strings of sites extending to the site j = 1. We apply the Wick-Bloch-de Dominicis theorem 
and present the result as the Pfaffian of the 2(2 j+n — 1) x 2(2j + n — 1) antisymmetric matrix 
constructed from the known elementary contractions (only these quantities are influenced by 
the existence of three-site interactions). Finally, we evaluate numerically the Pfaffians ob- 
taining as a result the desired xx time-dependent spin correlation function. To get S xx (k, to) 
(jl.8p we perform numerically the integration over time t and then the summation over n. 
Typically we take N = 400, assume j = 41, calculate (s|(t)s| +n ) for n up to n max = 100 
in the time range up to t max = 100. (However, for large fl we assume j = 81, n max = 100, 
^max = 200 whereas for (3 = 0.1 it is sufficient to take j = 41, n max = 50, t max = 50.) In 
the low-temperature strong-field regime the xx time-dependent spin correlation functions 
display long-time oscillations which lead to evident problems with integrating over time; 
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FIG. 10: S xx (k,uj) for the model <H3} with J= 1, D = E = 0, K = 0.5, O = (panel a), K = 2, 
n = (panel b), = 2.5, = (panel c), if = 0.5, O = 1.25 (panel d), K = 2, O = -1.125 
(panel e), K = 0.5, = —0.75 (panel f), K = 0.5, = —1.125 (for this set of parameters we 
assume e = 0.02, see the main text) (panel g) at the low temperature [3 = 20. The panels h and i 
correspond to high-temperature limit, (3 = 0.1, K = 0.5 and O = (h) or O = —1.125 (i). 



therefore, in this case (in fact, only for the set of parameters corresponding to panel g in 
Fig. [10]) we introduce under the integral in Eq. (11. 8p an auxiliary damping factor exp(— e|i|), 
where e is a small positive number. We examine in detail different types of finite size effects 
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401 ] to be sure that our results for S xx (k,u) refer to the thermodynamic limit. In Fig. 



[TU1 we demonstrate the results for the xx dynamic structure factors for a set of parameters 
which is in correspondence with the one used in Fig. [2] (panels a - f) and the analytical 
predictions (14.61) . (14.71) (panels g, h, i). 

We start to discuss the results obtained for S xx (k,uj) from the case of low temperatures. 
Although S xx (k, uj) is a many-fermion quantity and therefore is not restricted to a certain 
region in the k - uj plane it is mostly concentrated along certain lines in the k - uj plane 
which roughly correspond to the characteristic lines of the two-fermion excitation continuum 
discussed in Sec. ITT] (compare Figs. HUJi, [TJJb, [10b with Figs. EK [2)o, [2b (and Figs. EK EH, 
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Eg) as well as Figs. (TUft, US?, [TUf with Figs. EH, EK EF). For example, for J = 1, K = 0.5, 
f2 = (panel a in Fig. [TO]) S xx (k, uj) is accumulated along the three lines in the k - uj 
plane, f2_(/t) (13. 4p . uj^(k) 03.81) . shifted along K-axis by tc. Although the important role of 
the two-fermion excitations for the low-temperature many-fermion dynamic quantities like 
S xx (k,u) was noted several times earlier, we still do not have a simple explanation for that 
fact. On the other hand, the panels d, e, f and g in Fig. [10] demonstrate the development 
towards the zero-temperature strong- field result (14 .7p . 

As temperature increases the role of the two-fermion excitations diminishes and in the 
high-temperature limit many-fermion excitations produce the ^-independent Gaussian decay 
(compare Figs, [ffih, EE to Figs. QUi, [JUg). 



V. CONCLUSIONS 



In conclusion, we have examined several dynamic structure factors of the spin-1/2 trans- 
verse XX chain with (XZX + Y ZY)-tjpe three-site interactions. These three-site interac- 
tions essentially enrich the ground-state phase diagram of the spin model which may show 
two different spin liquid phases (spin liquid I and spin liquid II) in addition to the fer- 
romagnetic phase. We have calculated explicitly several dynamic structure factors (with 
the transverse dynamic structure factor S zz (k,uj) and the dimer dynamic structure factor 
Sj(k,uj) among them) which are governed exclusively by two-fermion (particle-hole) ex- 
citations. We have discussed in some detail the properties of the two-fermion excitation 
continuum determining its boundaries, soft modes and exponents of van Hove singularities 



4=11 ] . We have also discussed some specific features of different two-fermion dynamic struc- 



ture factors. Our analysis of many-fermion dynamic structure factors is restricted to the xx 
dynamic structure factor S xx (k, uj). For this quantity we have reported exact analytical re- 
sults in the high-temperature and zero-temperature strong-field limits and precise numerical 
results for other sets of parameters. The three-site interactions introduced leave a number of 
signatures in the dynamic quantities producing an extra two-fermion excitation continuum, 
van Hove singularity with exponent 2/3, or singularity of the dimer structure factor along 
the upper boundary of the two-fermion excitation continuum. In the presence of three- 
site interactions the symmetry of the xx dynamic structure factor in the zero-temperature 
strong-field regime with respect to the change Q — > — Q is broken. 
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We emphasize that the advantage of the model considered is its exact solvability, that 
means, in particular, the possibility to calculate various dynamic quantities accurately On 
the other hand, although there are some examples of real-life systems which can be modeled 
as spin-1/2 XX chains (see, e.g., Ref. 1421 ). the three-site interactions introduced are of a 
rather special kind, however, the reported results may serve to test other (approximate) 
techniques used to study more realistic models, e.g., with next-nearest neighbor interactions 
or with four-site interactions. Moreover, our results on dynamics may be used for discussing 
the effects of stationary energy fluxes in quantum spin chains. Thus in Ref. yc the transverse 
(zz) dynamic structure factor for a model with D = K = was discussed in relation to 
possible experimental observation of energy- current carrying states in quantum spin chain 
compounds. 
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APPENDIX: TIME-DEPENDENT xx SPIN CORRELATIONS IN THE ZERO- 
TEMPERATURE STRONG-FIELD REGIME 



We can rewrite Eq. ( l4~4l) (Q > [43]) as follows 

4W*K+»> = exp (-i (n + y) tj C n (A.l) 



with 



1 r w 

C n = — / d/t cos(ftn) exp (— it ( J cos k — i^cos 2 k)) . 
n Jo 



(A.2) 



The function C n introduced in (1A.2j) for even n can be expressed in terms of the function 
$3(/3, 7, x, y) (see section 5.7.1 in Ref. |44j) 



i jh 2 

C Q = exv(iKt)$ 3 [-,l;-iKt,— 



C 2m = exp(iift) f>l) w^ ^^Jl , (l + \,l + 1; -iKt, ^f) : (A.3) 



r(m-/ + i)r 2 (/ + i) 
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here T(n) is the gamma function. 
We notice that in the case K — 



C n = (-i) n J„(Jt), (A.4) 

whereas in the case J = 

c= { exp (if) i-J n (f ) , n = 2m, (a 
jo, n = 2m + 1; 

here J n (^) is the Bessel function of the first kind. 

Finally, we notice that the long-time asymptotics for (s*(t)s* +n ) in the zero-temperature 
strong-field regime can be calculated accurately using the stationary phase method 
\K/J\ < 1/2 we have 

t^oo 1 , , 

C n « — exp (lKt) 

y^\J\i ^ ' 

x | exp (-Ut + if sgn(J)) + exp (iJt - if sgn(J) + inn) ( ^ ^ ) 

For |tf/J| > 1/2 we have 



45|. For 



C n ta - expjiKt) 
^/2ir\J\t 

exp (-i Ji + if sgn( J - 2if)) exp (iJt - if sgn( J + 2K) + i?m) 
exp(-i(ir + 5)i+ifsg n (2if-i)) / j 



+2 = — — cos I narccos — - . (A. 7) 

1 2jc j I V 2K ' 1 



J 2K I 



For \K/J\ = 1/2 we have 



£ n * ~°° (_l)f( 1 -^)) ( F ^ l 6 xp (isgn(JK) (-\jt + ^sgn(J) 



2tt (2\J\t) : 



2 



' exp ( isgn(JK) ( -Jt - ^-sgn(J) ) +mn ) I . (A. 8) 



2^/nJ7\i V. V 2 4 

The long-time asymptotic behavior ( 1A.6I) . ( 1A.7I) . ( 1A.8I) may emerge already at relatively 
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short times as can be seen in Fig. [91 
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